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In the present work we focus on the case of (few-site) configurations respecting the PT-symmetry. 
We examine the case of such "oligomers" witli not only 2-sites, as in earlier works, but also the cases 
of 3- and 4-sites. While in the former case of recent experimental interest, the picture of existing 
stationary solutions and their stability is fairly straightforward, the latter cases reveal a considerable 
additional complexity of solutions, including ones that exist past the linear PT-breaking point in the 
case of the trimer, and more complex, even asymmetric solutions in the case of the quadrimer with 
nontrivial spectral and dynamical properties. Both the linear stability and the nonlinear dynamical 
properties of the obtained solutions are discussed. 



INTRODUCTION 

Over the past decade, the examination of Hamiltonian 
nonlinear dynamical lattices, as well as that of continuum 
systems with periodic potentials has been a subject of in- 
tense investigation The motivation for such studies 
stems from a variety of physical settings including, among 
others, the themes of optical beam dynamics in coupled 
waveguide arrays or optically induced photonic lattices 
in photorefractive crystals [2[ , the temporal evolution of 
Bose-Einstein condensates (BECs) in optical lattices Si , 
or the DNA double strand denaturation in biophysics J] . 
One of the common focal points among all of these ar- 
eas has been the intense study of the existence, stability 
and dynamical properties of their nonlinear (often local- 
ized in the form of solitary waves) solutions which are of 
principal interest and experimental observability within 
various applications; see [5[ for a relevant recent review. 

On the other hand, as the understanding of the conser- 
vative aspects of such systems comes to a point of mat- 
uration, a number of interesting variants thereof arise. 
A canonical one concerns the examination of effects of 
damping and driving that not onlyyield novel theoret- 
ical solutions (see as an example [6|), but also are in- 
herently relevant to applications (again, sec for a recent 
example Q])- A more exotic variant which, however, in 
the past couple of years has gained considerable momen- 
tum especially due to the recent experiments of Q is 
that of PT-symmetric dynamical lattices. This theme 
follows the pioneering realization of Bender and cowork- 
ers 0] that non-Hermitian Hamiltonians can still yield 
real spectra, provided that they respect the Parity (P) 
and Time-reversal (T) symmetries. Practically, in the 
presence of a (generally complex) potential the relevant 
transformations imply that the potential satisfies the con- 
dition V{x) = V*{—x). In nonlinear optics, the interest 
in such applications was initiated by the key contribu- 
tions of Christodoulides and co-workers [l^l which con- 
sidered solitary waves as well as linear (Floquet-Bloch) 
eigenmodes in periodic potentials satisfying the above 
condition, also including the effects of Kerr nonlinear- 



ity and observing how the properties of such waves were 
modified by genuinely complex, yet PT-symmetric po- 
tentials. 

More recently, motivated by the experimental possi- 
bilities and the relevant realization of a PT- "coupler" in 
0, there has been an interest in merging the experience 
of the above two areas, namely the consideration of PT- 
symmctric settings but for genuinely discrete media. In 
that vein, the experimentally-probed two-site system has 
been considered in the work of [llj, where it was shown 
that it can operate as a unidirectional optical valve, as 
well as in the study of [l^l, where the role of nonlin- 
earity in allowing (if sufficiently weak) or suppressing (if 
sufficiently strong) time reversals of exchanges of opti- 
cal power between the sites. Another recent example 
consisted of the generalization of where a lattice of 
coupled gain-loss dimers was considered. This theme has 
also been considered in the BEG literature and in the 
context of the so-called leaky Bose-Hubbard dimers (al- 
lowing e.g. the tunneling escape of atoms from one of the 
wells of a double- well potential). There, a variant of the 
model considered below has been self-consistently derived 
in the mean-field approximation [14| and the correspon- 
dence of its classical with the full quantum behavior has 
been explored [l^. 

Our aim in the present work is to revisit the examina- 
tion of the PT-symmetric coupler and to give a simple 
and complete characterization of the existence and sta- 
bility properties of its stationary solutions. It should be 
noted that this aspect has been partly addressed in both 
11 1 and [l2l- Nevertheless, we aim to give a characteri- 
zation thereof as a preamble towards the more complex 
(and thus, arguably, more interesting) generalization to 
what we call "PT-symmetric oligomers" , namely the con- 
sideration of a PT-symmetric trimer and that of a PT- 
symnictric quadrimer. Our aim here is to explore how 
the complexity of the problem expands as more sites are 
added, in order to offer a glimpse how such oligomers 
gradually give way to the elaborate phenomenology of a 
PT-symmetric lattice. We illustrate, for example, how it 
is possible in the case of a trimer to identify stationary so- 
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lutions which exist past the limit of hnear PT-symmetry 
breaking (something which is not possible in the dimer 
case) . We then proceed to illustrate how the phenomenol- 
ogy of the quadrimer is even richer and more complex, 
featuring among others asymmetric solutions with a re- 
duced symmetry spectrum differently than is the case for 
both the dimer and trimer. 

Our presentation will be structured as follows. In sec- 
tion II, we consider the fundamental (and previously con- 
sidered) dimer case. We use this as a benchmark for the 
presentation of our methods and results. Wc then turn to 
the more complex trimer case in section III and conclude 
our results with section IV on the quadrimer. Finally, 
section V summarizes our findings and presents some in- 
teresting questions for further study. 

DIMER 

We start our considerations from the so-called PT- 
symmetric coupler or dimer (as we will call it hereafter). 
In this case, the dynamical equations arc of the form: 

iui = —ku2 — juipui — i'-fui 

iu2 = — fcui — |u2pU2 + i7U2. (1) 

The model of Eq. ([1]) considers the linear PT-symmetric 
dimer experimentally examined in Q, as augmented by 
the Kerr nonlinearity relevant e.g. to optical waveguides; 
see also [III ■ The overdot denotes the derivative with 
respect to the evolution variable which in optical appli- 
cations is the propagation distance. In what follows, we 
will denote this variable by t (to indicate its evolution- 
ary nature). We seek stationary solutions of the form 
Ml = cxjp{iEt)a and U2 = exjp{iEt)b. Then the station- 
ary equations arise: 

Ea ~ kb + \a\'^a + i^a 

Eh = ka+ \h\H~i-fh. (2) 

Using a generic polar representation of the two "sites" 
a = Ae'^'^" , b = Be^'^'' , we are led to the following alge- 
braic conditions for the two existing branches of solutions 
(notice the ± sign distinguishing between them): 

^B"^ ^E± v/fc2 - 72 (3) 
sin((/)b- 0„) = (4) 

The fundamental difference of such solutions from their 
standard Hamiltonian (7 = 0) counterpart is that the 
latter were lacking the "flux condition" of Eq. (|4]). This 
dictated a selection of the phases so that no phase current 
would arise between the sites. On the contrary, in PT- 
symmetric settings, the phase flux is nontrivial and must, 
in fact, be consonant with the gain-loss pattern of the 
coupler. 
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FIG. 1: The two branches of solutions for the dimer problem 
are shown for parameter values k — E — 1. The top left illus- 
trates the amplitude of the sites, the top right their relative 
phase and the bottom panel shows the (nontrivial) squared 
eigenvalue of the two branches. The solid line corresponds to 
the always stable branch u'^' , while the dashed branch corre- 
sponds to it'^\ which acquires a real eigenvalue pair above a 
certain 7 = — E'^/A. 



Fig. [T] shows the proflle of the two branches. The flrst 
branch u^^^ corresponding to the (— ) sign in Eq. ^ 
is stable when 7^ < fc^ — E'^/A, whereas the second 
branch u^^) is always stable. The linearization around 
these branches can be performed explicitly yielding the 

nonzero eigenvalue pairs ±2i-^2(fc2 — 7^) — E^/k^ — 7^ 

for the first and ±2i^2{k^ ^ 7^) + E\fk'^ — 7^ for the 
second (notice that the latter can never become real). 

It is relevant to note here that the two branches "die" 
in a saddle-center bifurcation at 7 = fc, as shown in the 
figure. Importantly, this coincides with the linear limit 
of the PT-symmetry breaking since the linear eigenvalues 
of the problem are A = i-^/fc^ — 7^. Hence, the nonlin- 
ear solutions terminate where the linear problem eigen- 
functions yield an imaginary pair, predisposing us for an 
asymmetric evolution past this critical point (for all ini- 
tial data) . The dynamical evolution of the dimer is shown 
first for a case of 7 < (in which is unstable, while 
is stable) in Fig. [5] The evolution of the instability 
of u^^^ leads to an asymmetric distribution of the power 
in the coupler, despite the fact that parametrically we 
are below the linear critical point (for the PT-symmetry 
breaking). Notice that in all the cases, also below, where 
a stationary solution exists for the parameter values for 
which it is initialized, dynamical instabilities arise only 
through the amplification of roundoff errors i.e., a numer- 
ically exact solution up to 10 ~® is typically used as an 
initial condition in the system. Naturally, beyond 7 = fc, 
as shown in Fig. |31 all initial data yield such an asym- 
metric evolution. 
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FIG. 2: Dynamical evolution of initial data belonging to the 
two branches of stationary solutions of a dimer in the case of 
7 = 0.9, E = k = 1, which is past the critical point for the 
instability of the first branch (left panel), while the second 
branch of the right panel is still dynamically stable. Notice 
that the left panel is plotted in semilog. 
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FIG. 3: Similar initialization of the dimer based on the two 
branches of stationary solutions (for 7 = 1) but now for the 
case of 7 = 1.1 {E = k = 1). The asymmetric evolution of 
the coupler past the linear PT-symmetric threshold can be 
clearly discerned. 



TRIMER 

We now turn to the case of the trimer where the dy- 
namical equations arc 



lUl 
iU2 

ills 



-ku2 — \ui\^Ul — l^Ux 
-k{ui + U3) - |w2pU2 



-ku9 



(5) 



Seeking once again stationary solutions leads to the al- 
gebraic equations 



Ea = kb + \a\ a + ija 
Eh ^ k{a + c) -f \h\'^h 
Ec ~ kb + |cpc — ijc 



(6) 



In this case too, it is helpful to use the polar repre- 
sentation for the three-sites in the form a = ^e*"^" , b = 
Be^'^'',c ~ Ce^'^", which, in turn, leads to the algebraic 
equations of the form: 



sm 



B'^ - EB^ + 2EA^ -2A^ = Q 

jA 

h - 4>a) = - sin(0b - (j)c) = 

EA-A^ 



cosi 



cos 



kB 



(7) 

(8) 

(9) 
(10) 



Notice how the presence of the gain-loss spatial profile 
along the 3-sites induces a spatial phase distribution and 
enforces the condition of a symmetric amplitude profile 
with the first and third site sharing the same amplitude. 
This phase distribution would be trivial (relative phases 
of or tt) in the 7 = case. 
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FIG. 4: Existence and stability of solutions for the case of the 
trimer analogously to Fig.[T]for the dimer case. The only other 
difference is the illustration of both real and imaginary parts 
for the eigenvalues in the bottom panel. We have set (without 
loss of generality) c^i, normalized to 0. Different line styles are 
used to distinguish different branches and their corresponding 
stabilities (see the detailed explanation in the text). 

A typical example of the branches that may arise in 
the case of the trimer is shown in Fig. |4]for = fc = 1. 
In this case, we find three distinct branches in the consid- 
ered interval of parameter values. There are two branches 
which exist up to the critical point 7 = 1.043. In this in- 
terval one of the two branches u^^^ is chiefly unstable 
(denoted by dash-dotted line) except for a small interval 
of 7 e [1, 1.035]. The other one u*-^^ is chiefly stable (de- 
noted by a dashed line) except for 7 e [1.035, 1.043]. The 
eigenvalues of u^^'' and u*^^^ in 7 S [1.035, 1.043] are very 
close to each other but not identical. Notice that u^^"^ 
is unstable due to a complex eigenvalue quartet whose 
eigenvalues collide on the imaginary axis for 7 = 1 and 
split into two imaginary pairs one of which becomes real 
for 7 > 1.035. Finally, these two branches collide in a 
saddle-center bifurcation (for 7 = 1.043) and disappear 
thereafter. 

Interestingly, however, these are not the only branches 
that arise in the trimer case. In particular, as can be seen 
in Fig. m there is a branch of solutions bifurcating from 
zero (amplitude) for 7 > y/2k'^ — E'^, denoted by u'^^\ the 
solid line in Fig. |4l In our case E = k = 1, this branch is 
only stable for 7 < 1.13, at which point two pairs of imag- 
inary eigenvalues collide and lead to a complex quartet 
which renders the branch unstable thereafter. Yet, this 
branch of solutions has a remarkable trait. In the case of 
the trimer, the underlying linear problem possesses the 
following eigenvalues 0, ±1/2^^ — 7^. Hence, the criti- 
cal point for the existence of real eigenvalues of the lin- 
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ear problem in the case of the PT-symmetric trimer is 
7 = \f2k (cf . with the 7 = A; Hmit of the dimer) . Never- 
theless, and contrary to what is the case for the dimer, 
the third branch of solutions considered above persists 
beyond this critical point (although it is unstable in that 
regime). 

The evolution of the three distinct branches of solu- 
tions, namely the chiefly unstable one the chiefly 
stable one u^^^ and flnally of the one persisting past the 
linearly unstable limit u^^^ is shown, respectively, in Figs. 
[51 ini and [71 It can be seen that in accordance with 
the predictions of our linear stability analysis the first 
two branches are stable or unstable in their correspond- 
ing regimes, while past the point of existence of these 
branches (7 = 1.043) their evolution gives rise to asym- 
metric dynamics favoring the growth of the power in a 
single site (or in some cases even in two sites; see e.g. 
the bottom panels of Figs. [5]and|6l). On the other hand, 
for the branch emerging at 7 = 1 and persisting past 
the linear instability limit, we indeed find it to be stable 
for 7 < 1.13 and unstable thereafter again leading to an 
asymmetric distribution of the power. 
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FIG. 6: Similar to Fig. [5] but for the case of vP'^ and for 
E = k = 1. The bottom panel is again initialized for 7 — 1.04. 
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FIG. 7: The dynamical evolution for the case of u^^^ and for 
E — k = 1; ^ = 1.04 is used in the left panel and 7 = 1.5 in 
the right one. 
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FIG. 5: The dynamical evolution of the amplitudes of the 
three sites of the stationary solution u*-^^ in the case oi E = 
k = 1 for 7 = 0.5 (left panel), 7 — 1.03 (right panel) and 
7 = 1.1 (bottom panel). The bottom panel is initialized with 
the exact stationary solution for 7 = 1.04 (since for 7 = 1.1 
the branch no longer exists as a stationary solution). 



Notice here that we only consider the case where the first 
two sites have the same loss and the latter two the same 
gain. This is by no means necessary and the gain-loss 
profile can be generalized to involve two-parameters (e.g. 
±71 and ±72 distinct between the different correspond- 
ing sites i.e., the first and fourth ones, as well as the 
second and third ones). We do not consider this latter 
case here, due to its more complicated algebraic structure 
that does not permit the direct analytical results given 
below. More specifically, in our considered special case, 
the stationary equations read: 



QUADRIMER 

Finally, we briefly turn to the case of the quadrimer. 
Here the equations are: 

iui = — fcM2 — Iwipui — 47^1 

iU2 = -fc(Ml -t- U3) - |u2pM2 - *7W2 

ills = -k{u2 + U4) ~ \u3\'^U3 + ijU3 

iU4 = —ku3—\u4\'^U4 + i-fU4 (11) 



Ea = kb + \a\'^a + i'ja 

Eh = fc(a + c) + \b\% + i-ih 

Ec = k{b + d) + \c\^c-ijc 

Ed = kc+\d\^d~ijd (12) 

The polar representation of the form a = Ae'^'^" , b = 
Be^^^c = Ce*'^'=,d = De^'^'' now allows the following 
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reduced algebraic equations: 



= E 

A^B^ + j^A^ - k^B^ = 



sm 



sin 



b) = - 



sm 



kBC 



= -1 

kC' 



(13) 
(14) 
(15) 

(16) 
(17) 
(18) 



Notice that in this case not only do we have the custom- 
ary phase profile, but in fact one of the phase differences 
becomes locked to 7r/2 due to the presence of the gain- 
loss pattern. 

Upon reducing the algebraic equations, we obtain 

[E - B^)B'^ +j^{E - B^) - k^B^ ^0 (19) 
{E ^C^)C'^ + j^{E -C^} - k^C^ ^0 (20) 

jE = kBC (21) 

This leads to the important conclusion that for this gain- 
loss profile in the case of the quadrimer, differently than 
in the cases of the dimer and trimer, one of the parame- 
ters E, fc, 7 is determined by the other two; i.e., not aU 
three of these parameters can be picked independently in 
order to give rise to a solution of the quadrimer. 

We hereby set E ^ 1, and increase 7 from as before, 
then k can be obtained self-consistently from the above 
equations. Therefore, once E and 7 are fixed, the solu- 
tions of the quadrimer problem are fully determined. We 
now present three branches of solutions that arise in this 
setting, as we increase 7. These are shown in the panels of 
Fig. [8] There are two classes of solutions here. The solid 
curve u'^' corresponds to a fully asymmetric branch with 
A, B, C, D distinct, something that is unique (among 
the settings considered herein) to the quadrimer. Fur- 
thermore, this always unstable branch does not respect 
the Hamiltonian eigenalue symmetry i.e., that if A is an 
eigenvalue, so are —A, A* and —A*. On the other hand. 



the dashed curve of the branch 



and the dash-dotted 



curve of u"' correspond to symmetric branches with am- 
plitudes A = D and B — C . Among the two symmetric 
branches m^^-* and m^'^-' that collide and disappear together 
in a saddle-center bifurcation at 7 = 0.362, we can ob- 
serve that the former one between them has a real and 
two imaginary pairs of eigenvalues being always unsta- 
ble, while the latter starts out stable, but the collision of 
two of its imaginary pairs will render it unstable past the 
critical point of 7 = 0.023. Interestingly the asymmet- 
ric branch u^^^ and the symmetric branch u^^'^ appear to 
collide in a subcritical pitchfork bifurcation that imparts 
the instability of the asymmetric branch to the symmet- 
ric one for 7 > 0.193. 



As an aside, we should also note here that in its linear 
dynamics (examined e.g. experimentally in Q) the PT- 
symmetric quadrimer has an interesting difference from 
the dimer and trimer. In particular, the 4 linear eigen- 
values of the system are: 



Ai,2=±y-7' + -(3fc-y^T6^2^5fc2) (22) 

A3,4 = + ^ (^^ + + 5fc2) . (23) 

The fundamental difference of this case from the others 
considered above is that these eigenvalues do not become 
imaginary by crossing through 0. Instead, they become 
genuinely complex, through their collision which occurs 
for 7 = •\/5fc/4, a critical point which is lower than that 
of the trimer. This could be an experimentally observ- 
able signature of the difference between the near linear 
dynamics of the quadrimer in comparison e.g. to the 
trimer. 
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FIG. 8: Three branches of sohitions for the trimer problem 
with parameters E = 1 and (j)a normalized to 0: the solid 
lines denote the asymmetric branch m'^', while the dashed 
and dash-dotted denote the symmetric branches u^'^^ and w'-'^' , 
respectively. 

The dynamics of these different branches was also con- 
sidered in Fig. ini In this case, it can be clearly observed 
that all three branches tend towards an asymmetric dis- 
tribution of the power. This favors the two sites (third 
and fourth) with the gain, although some case examples 
can be found (see e.g. the top left panel of Fig. [S] for the 
asymmetric branch) , where only one of the two gain sites 
is favored by the mass evolution. 



CONCLUSIONS 

In the present work, we considered the existence, sta- 
bility and dynamics of PT-symmetric oligomers i.e., con- 
figurations with few sites. Similarly to the recent works 
of 
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12l | and also the experimental investigation of [8|] , 
we have started our considerations by a complete char- 
acterization of the dimer case, where the two obtained 
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FIG. 9: The profile of the dynamical evolution of the three 
different branches (top left for u^^\ top right for u'^^ and 
bottom panel for u'^') of a quadrimer in the case of i5 = 1 
and 7 — 0.1. 

branches of solutions terminate at the critical point of 
the linear case. However, we illustrated that the trimer 
and quadrimer feature a number of fundamental differ- 
ences in comparison with this dimer behavior. In partic- 
ular, the trimer features branches which exist past the 
linear critical point (although unstable). On the other 
hand, the quadrimer has even richer features: in partic- 
ular, it possesses asymmetric solutions whose spectrum 
only has symmetry around the x-axis (and not the four- 
fold symmetry of the 7 = Hamiltonian problem) . The 
bifurcation structure is also richer in the latter problem 
featuring symmetry-breaking pitchfork bifurcations. An- 
other notable feature is that solutions do not exist for 
arbitrary combinations of coupling, gain/loss parameter 
and propagation constant; instead, these parameters ap- 
pear to be inter-connected (at least in the case of a single 
gain- loss parameter considered herein). Finally, even the 
linear problem presents interesting variations in this case, 
featuring the breaking of the real nature of the eigenval- 
ues through two colliding pairs that lead to a quartet 
occuring for smaller gain/loss parameter values than in 
the trimer case. 

This investigation may be a first step towards obtain- 
ing a deeper analytical understanding of the features of 
PT-symmetric lattices. In such settings it would be rel- 
evant to obtain general conclusions both for the linear 
dynamics (and how it depends on the gain/loss profile 
parameters) as well as more importantly for the nonlin- 
ear modes, including the solitary waves that may arise. 



Understanding such modes and the comparison of their 
properties to the continuum ones, as well as to the dis- 
crete ones in the absence of the gain/loss would be im- 
portant directions for future study. 
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